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We study evolution of quark-gluon matter in the ultrarelativistic heavy-ion collisions within the 
frame work of relativistic second-order viscous hydrodynamics. In particular, by using the various 
prescriptions of a temperature-dependent shear viscosity to the entropy ratio, we show that the 
hydrodynamic description of the relativistic fluid becomes invalid due to the phenomenon of cavi- 
tation. For most of the initial conditions relevant for LHC, the cavitation sets in very early stages. 
The cavitation in this case is entirely driven by the large values of shear viscosity. Moreover we also 
demonstrate that the conformal terms used in equations of the relativistic dissipative hydrodynamic 
can influence the cavitation time. 
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Presently the viscosity of the strongly-interacting mat- 
ter produced in the heavy-ion collision experiments at 
LHC and RHIC is under extensive investigations. The 
measurements of the elliptical flow parameter u 2 show a 
strong collectivity in the fluid-flow implying the existence 
of a very low viscous stress due to shear viscosity pQ . Ac- 
cording to the AdS/CFT conjecture, ratio of the shear 
viscosity to entropy density rj/s may not be lower than 
1/47T which is now known as KSS bound [2J. It has been 
argued that in order to explain the collective flow data 
•q/s cannot be larger than twice the KSS-bound [3]. 

It must be noted that the applications of the viscous 
hydrodynamics discussed above regard rj/s as indepen- 
dent of temperature. However, recently it has been ar- 
gued that constant rj/s is in sharp contrast with the 
observed fluid behavior in nature where it can depend 
on temperature [HE]. It has been demonstrated that 
the temperature-dependence of r\j s can strongly influ- 
ence the transverse momentum spectra and elliptical flow 
in the heavy-ion collision experiments at LHC [SJ [5] . It 
should be emphasized here that the ratio of bulk viscos- 
ity to entropy density £/s as a function of temperature 
was already considered by several authors and interest- 
ing consequences like cavitation were studied [7HTU], A 
similar analysis with a temperature-dependent rj/s has 
not been performed so far, which we intend to address 
here. Cavitation has also been studied recently with a 
holographic formulation of sQGP [IT] . 

It is generally expected that 77/s for QGP has a min- 
imum at the critical temperature T c , while it increases 
with the temperature beyond T c Q21 [13]. In this 
work we use rj/s prescriptions arising from lattice QCD 
(1QCD) as in Ref. 5j, virial theorem type of arguments 
[12) as well as the analytical expressions for 77/s as given 
in Ref. jB]. We show that the large values of 77/s, relevant 
for LHC energies, can make the effective pressure of the 
fluid very small in a time less than 2 fm/c. This would 
cause cavitation in the fluid which in turn would limit 
the applicability of hydrodynamics. It must be noted 
that the cavitation at RHIC energies studied in Refs. 
[SI HI [14] earlier was driven by the high values of the 
bulk viscosity near the critical temperature. However, 
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FIG. 1. Different prescriptions of rj/s as function of temper- 
ature, with T c =0.2 GeV. The horizontal curve show r)/s = 
1/471" obtained from the AdS/CFT correspondence. 



the bulk viscosity can play an insignificant role in the 
temperatures T » T c . In the present study we demon- 
strate that for LHC energies cavitation is solely driven 
by the shear viscosity. 

We use relativistic boost invariant causal viscous 
hydrodynamics equations in 1+1 dimensions [151 116] . 
One may argue against the validity of applying (1+1)- 
dimensional flow in studying the relativistic heavy-ion 
collisions by ignoring the transverse flow. As will be 
shown later for a central collision at LHC energies the 
cavitation sets during the initial stage of the evolu- 
tion in a time less than 2 fm/c. Since the transverse 
flow is negligible during the earlier stages of a heavy- 
ion collision, it will not have a significant effect on the 
cavitation time. We use the parametrization of the 
coordinates t = r cosh?7 s and z = r sinh?7 s , with 
the proper time r = \/t 2 — z 2 and space-time rapid- 
ity j] s = 5^i[xzJ]. Now the 4-velocity can be written 
as u M = (cosh?7 s , 0, 0, sinh?7 s ). Within the second order 
theory (for more details on this theory and its applica- 
tion to relativistic heavy ion collisions we refer to Refs. 
[S][T7][IH]) the equations dictating the longitudinal ex- 



pansion of the medium are given by 
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The effects due to shear and bulk viscosity are repre- 
sented via <E> and II respectively and they can contribute 
to the effective pressure of the fluid. Eqns. (2][3) are evo- 
lution equations for $ and II governed by their relaxation 
times tv and Tp respectively. Last term on the right-hand 
side of Eq. ^ is due the conformal symmetry [53]. In 
order to close the system of Eqns. (l][3), one needs to 
use equation of state (EoS). We have used recent 1QCD 
results [24] for this purpose. At LHC energies the bulk 
viscosity is expected to be negligible as e ~ 3P and one 
can ignore Eq. (§. 

In the local rest frame the shear stress describes the de- 
viation from the isotropy of the stress tensor. To quantify 
this anisotropy one can define the longitudinal pressure 
P z [HI [23 in absence of the bulk stress as 



P, = P - $ 



(4) 



where P is the equilibrium hydrodynamic pressure. 

We use recent 1QCD estimate for rj/s in QGP sector 
calculated by Nakamura et. al |13| . The resulting rj/s 
from 1QCD has the expected minimum near the criti- 
cal temperature T c . It should be noted that recent lat- 
tice studies indicate a crossover rather than a phase- 
transition |26| . However, for the present work this may 
not be an issue since we are interested in temperature 
dependence of rj/s where T c is a parameter. We use 
the parametrization of rj/s given in Ref. [37], where 
the minimum value of rj/s is I /An. Another prescrip- 
tion for shear viscosity that we use is from Ref. [12], 
where using virial expansion techniques, the authors cal- 
culate rj/s in QGP. Fig.[l] shows the plots of various 
rj/s prescriptions versus temperature with T c =0.2 GeV. 
The top curve shows values of rj/s obtained from the 
lattice results, while the middle curve corresponds to 
•q/s values obtained from the virial expansion. The 
horizontal line corresponds to the KSS value. Finally, 
we consider the temperature-dependent forms of rj/s as 



given in Ref. [BJ: 
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TU65 GeV. 

Relaxation time tv — 2t7/3 2 can be determined by an 
underlying theory other than the hydrodynamics. It 
ought to be mentioned that in the relativistic viscous hy- 
drodynamic literature there is some ambiguity regarding 
the value of the relaxation times associated with shear 
and bulk viscosity. In this work we have taken the relax- 
ation time for shear viscosity t v = which is moti- 
vated by kinetic theory [SI [57J ■ In addition we also solve 
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To=0.3 fm/c 
T =0.506 GeV 


rj/s 1QCD 


21.18 


0.57 


0.421 


14.21 


0.52 


0.434 


rj/s virial 


20.61 


0.63 


0.410 


13.93 


0.93 


0.374 


to=0.3 fm/c 
T =0.560 GeV 


rj/s 1QCD 


31.58 


0.57 


0.465 


20.31 


0.52 


0.479 


rj/s virial 


25.06 


0.68 


0.444 


18.12 


1.20 


0.385 


to=0.6 fm/c 
T =0.405 GeV 


rj/s 1QCD 


12.40 


1.20 


0.333 


10.83 


1.53 


0.316 


rj/s virial 


15.30 


1.27 


0.329 


11.88 






to=0.6 fm/c 
T =0.450 GeV 


rj/s 1QCD 


18.36 


1.21 


0.369 


15.63 


1.29 


0.365 


rj/s virial 


19.84 


1.43 


0.353 


16.07 






to = 1.0 fm/c 
T =0.350 GeV 


■q/s 1QCD 


10.48 






9.98 






r//s virial 


13.04 


2.16 


0.283 


11.17 







TABLE I. Column IS corresponds to the case when the con- 
formal terms are neglected from the hydrodynamics equa- 
tions. In this case the relaxation time tv from the kinetic the- 
ory is taken in to account. The column IS+C corresponds to 
the case when the conformal terms and tv obtained from the 
supersymmetric Yang-Mills theory are included in the equa- 
tions of hydrodynamics. The cavitation time r cav and t/ are 
measured in the unit of fm/c and the cavitation temperature 
T cav is shown in the units of GeV. T cav and T cav are left blank 
when there is no cavitation. 
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Eqns. (l][2) by taking tv — , , _ — , 

Ai = or' ms pi re d by results from Af=4 supersymmetric 
Yang-Mills theory [2^1135] , 

Next we present the numerical solutions for the equa- 
tions of hydrodynamics. First we consider the case with 
temperature-dependent rj/s taken from 1QCD calcula- 
tions. Fig. [2] shows the plots of longitudinal pressure P z 
versus the proper time for the cases of pure Israel-Stewart 
type (IS) hydro by neglecting the conformal terms in Eq. 
([2]) and with conformal terms (IS+C). In the case of IS 
we use tv from the kinetic theory and from the super- 
symmetric Yang-Mills theory when we consider IS+C 
case. We plot P z for these two cases with the initial 
temperatures 0.405 and 0.450 GeV. The starting time 
To is chosen to be 0.6 fm/c. Let us first consider the 
case with To = 0.405 GeV. From the figure it is clear 
that longitudinal pressure becomes negative in the IS 
case around cavitation time T cav = 1.20 fm/c. The tem- 
perature T cav at which the cavitation occurs is about 
0.333 GeV which is much larger than the critical tem- 
perature T c . Thus the cavitation can take place very 
early during the evolution. This, we believe, provides a 
posteriori justification for neglecting the transverse flow; 
as the hydrodynamic treatment may not be valid for the 
time larger than r cav . Further, if we include the con- 
formal terms in Eq. ^ together with the relaxation 
time obtained from supersymmetric Yang-Mills (IS+C), 
the cavitation time increases marginally and becomes 
T cau=L53 fm/c. Similarly T CQU =0.316 GeV is less than 
the cavitation temperature without the conformal terms. 
Next we consider a higher initial temperature To = 0.450 
GeV. Here also we observe cavitation for both IS and 
IS+C cases as in the previous case with T = 0.405 GeV. 
For IS case cavitation happens at a time r cav = 1.21 
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FIG. 2. The longitudinal pressure P z as function of time for 
IS and IS+C hydrodynamics. Initial time is taken to be 0.6 
fm/c with initial temperatures 0.405 and 0.450 GeV. r)/s(T) 
is obtained from the 1QCD curve shown in Fig.[l]. 



FIG. 3. The longitudinal pressure P z as function of time for 
IS and IS+C hydrodynamics. Initial time is taken to be 0.6 
fm/c with initial temperatures 0.405 and 0.450 GeV. rj/ s(T) 
taken from the virial expansion techniques curve in Fig. [I]. 



fm/c which is only marginally greater than the corre- 
sponding T = 0.405 GeV case considered previously. 
However, here the temperature at which cavitation oc- 
curs is higher (T cav = 0.369 GeV) than the previous 
case. This difference is expected since the initial tem- 
perature for the latter case is also larger. IS+C case 
with T = 0.450 GeV, cavitation sets in at r cav = 1.29 
fm/c with T cav =0.36b GeV. Again we note that there is 
not much difference between the cavitation times in IS 
and IS+C cases. 

In Fig. [3], we show P z as function of time by taking 
77/s values using the virial expansion techniques given in 
Ref. |12| . Values for r and T are same as in Fig. [2]. 
Here in the IS case with To = 0.450 GeV we can see 
that cavitation sets in around 1.43 fm/c when the system 
temperature is 0.353 GeV. However, as one can see from 
Fig. [3], when we include conformal terms (IS+C case) 
cavitation scenario is avoided. Next we lower the initial 
temperature to 0.405 GeV and consider the IS case. Here 
system reaches a negative longitudinal pressure stage at 
r cav = 1.27 fm/c with T cav = 0.329 GeV. But with con- 
formal terms included, as one can see from the figure, 
the longitudinal pressure remains positive although it as- 
sumes a very small value by 2 fm/c. Since the values 
of 77/s for the virial expansion techniques are systemati- 
cally smaller than 77/s for the 1QCD results as shown in 
Fig.yy, the corresponding cavitation time is larger than 
that shown in Fig. J2J. However, the cavitation temper- 
ature T cav is smaller than the corresponding cases dis- 
cussed in Fig. [2]. 

Further, we have changed the values of the initial 
time by considering the case To = 0.3 fm/c and To = 
1.0 fm/c. These results are summarized in Table I. For 
To = 0.3 fm/c and To = 0.560 GeV case, the cavitation 
occurs around r cav = 0.6 fm/c for the 1QCD 77/s while 
it occurs around r cav = 0.68 fm/c for 77/s obtained from 
the virial expansion. For the case with t = 1.0 fm/c 
and Tq = 0.350 GeV, for 77/s from virial expansion, the 



cavitation occurs around r cav — 2.16 fm/c. However, in 
this case when the 77/s values from 1QCD are used there 
is no cavitation. We would like to note that the table 
shows no entries for r cav and T cav for certain cases. For 
such instances the longitudinal pressure remains positive 
and there is no cavitation. Table I indicates for the given 
initial conditions there are more number of no-cavitation 
instances when the conformal terms in the equations of 
the hydrodynamics are taken into account. 

We also summarise the results for tj, the total time 
taken by the system to reach T c by ignoring the cavitation 
in Table I. One can see that with T = 0.405 GeV for 
1QCD (virial) case r/=12.40 (15.30) fm/c without the 
conformal term and t/=10.83 (11.88) fm/c if the term 
is included. Thus the inclusion of the conformal terms 
reduces Tf . We would like to emphasize that in this work 
we have taken a rather conservative initial value $(to)=0 
so that the initial value of the longitudinal pressure is 
always positive [29]. Instead if one includes the first- 
order (Navier-Stokes) initial value $(t ) = 4t7(T )/(3to), 
then the cavitation can occur at even earlier time and 
higher temperature. 

Next, we repeat our analysis using the temperature- 
dependent 77/s prescriptions given in Ref. [BJ. With the 
same initial conditions as in Ref. [6] we find that the 
longitudinal pressure becomes negative very early ~ 1 
fm/c for all the cases they have considered. Fig. [4] shows 
P z versus r for initial temperature To = 0.419 GeV and 
To = 0.6 fm/c. In this case also cavitation sets in early 
in about r cav ~ 1 fm/c. 

Further, we also consider 77/s = 1/4tt and £/s as func- 
tion of T as in Ref. [8] for LHC energies. It is found 
that the cavitation does not occur in this case unlike the 
results for RHIC energies [5] GO 03] ■ One may naively 
expect that when the system temperature reach T ~ T c , 
the bulk viscosity become large enough to drive cavita- 
tion. However, the cavitation occurs when the viscous 
stress (n and/or $) has a peak in its temporal profile 
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FIG. 4. Cavitation with various r]/s prescriptions considered 
by Shen et.al. in Ref. [6]. The initial temperature is taken to 
be 0.419 GeV with initial time 0.6 fm/c. 




FIG. 5. Cavitation along with anomalous viscosity. The lon- 
gitudinal pressure P z and <& as function of time. The initial 
temperature is taken to be 0.450 GeV with initial time 0.6 
fm/c. 



(see Fig. [5]). The height of the peak is determined by t , 
T and the initial values of C/s or r]/s. For LHC ener- 
gies, we find that even at the peak value of the viscous 
stress IT, the condition IT < P is satisfied and therefore 
cavitation does not occur. 

We have further considered the effect of anomalous vis- 
cosity (t/a), which may be important during the early 
time evolution in the hydrodynamics [30] . We use an ef- 
fective shear viscosity 77" 1 = r/^ 1 + 77^ 1 as discussed in 
Ref. [3D]. Here, r\c the collisional viscosity is taken from 
1QCD and for tja/s we use the expression from Ref. [3D]. 
In this case, (with To = 0.6 fm/c and Tq — 450 MeV), cav- 
itation sets in at a time 1.46 fm/c when the system is at a 
temperature 351 MeV. The initial value of anomalous vis- 
cosity to entropy density ratio is ~ 0.23. The results are 
presented in Fig. [5 , where we plot the shear stress term 
$ and longitudinal pressure P z as a function of proper 
time. As is clear from Fig. [5] the shear stress $ increases 
sharply from its initial value. The maximum value of <& 



and the time it takes to reach that value strongly depend 
upon . This sharp rise of $ result in a sharp reduction 
of P z , which, finally becomes negative at r cav . 

Thus to summarise, we have shown by using various 
prescriptions for a temperature-dependent rj/s that at 
LHC energies the higher values of shear stresses can in- 
duce the cavitation. This will in turn make the hydro- 
dynamic treatment invalid beyond cavitation time r cav . 
We have studied shear viscosity induced cavitation using 
one dimensional boost invariant causal dissipative hydro- 
dynamics of Israel-Stewart. One would of course like to 
do an analysis using a (3+l)-dimensional viscous hydro- 
dynamics like e.g. in Ref. [3T]. Since cavitation occurs 
during the early stages of the collision, we believe that 
the inclusion of transverse flow will not alter the result 
qualitatively. However, as a caveat, we would like to 
mention that the difference between the initial conditions 
for the "cavitation" and "no-cavitation" cases is rather 
small, see Table I. It remains to be seen if the inclu- 
sion of transverse flow can alter the cavitation scenario 
in a qualitative way. It is worth noting here the negative 
pressure scenario may be circumvented by considering 
anisotropic corrections in the distribution functions [35] . 
It should be emphasized that there exist alternate formu- 
lations of dissipative relativistic fluid dynamics where the 
longitudinal pressure remains positive e.g. in Ref. |33j 
It has been shown recently that the inclusion of the cav- 
itation condition in boost invariant hydrodynamics can 
change the particle spectra from expanding QGP [9|,ll4j. 
Based on the various prescriptions of r]/s our results indi- 
cate that the hydrodynamical description is valid about 
i~ C av ~ 2 fm/c at LHC energies. Beyond T cav , the fluid 
might fragment [7] or form inhomogeneous clusters. Let 
us note that one of the assumptions of the statistical 
hadronisation models lies in creation of extended clus- 
ters of quark matter which hadronize statistically |34) . 
Alternately, as has been attempted recently one can pos- 
sibly use a hybrid approach for the description of fire ball 
expansion applying viscous hydrodynamics for the QGP 
stage and then coupling it to a microscopic kinetic evo- 
lution for the hadronic stage [35]. Mere integration of 
the equations of hydrodynamics may not tell us about 
cavitation. We therefore believe that the conditions for 
cavitation may be required to be incorporated in the hy- 
drodynamical codes. 
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